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Abstract 

We provide a theory of dynamo (a effect) and momentum transport in three-dimensional mag- 
netohydrodynamics. For the first time, we show that the a effect is severely reduced by flow shear 
even in the absence of magnetic field. The a effect is further suppressed by magnetic fields well 
below equipartition (with the large-scale flow) with different scalings depending on the relative 
strength of shear and magnetic field. The turbulent viscosity is also found to be significantly re- 
duced by shear and magnetic fields, with positive value. These results highlight a crucial effect of 
shear and magnetic field on dynamo quenching and momentum transport reduction, with impor- 
tant implications for astrophysical and laboratory plasmas, in particular for the dynamics of the 
Sun. 

PACS numbers: 47.65.-d, 91.25.Cw, 95.30.Qd, 96.60.Hv 



1 



Introduction - Dynamo action describes the process through which a motion of a conducting 
fluid in the presence of a magnetic field amplifies that magnetic field. This is a fundamental 



mechanism that explains ubiquitous magnetic fie 



physical, geophysical and laboratory plasmas 



ds in a variety of systems, including astro- 



This is especially the case in stellar 



convective regions where the diffusive time-scale is so short due to turbulent motion that 
any primordial field would decay in a few million years. In a conducting fluid of velocity V, 
magnetic field B evolution is governed by the induction equation: 

d t B + V- VB = B- VV + i)V 2 B and V ■ B = , (1) 

where 77 is the ohmic diffusivity. In case of the Sun, the prominent radial shear layer 
(tachocline) permits the creation of a toroidal magnetic field from an existing poloidal one 
via shearing [the third term in Eq. (JT]), so-called Q effect]. The generation of poloidal 
field from the toroidal field requires another mechanism such as a effect via kinetic helicity, 
magnetic buoyancy [4] or some kind of magnetic instability [5] and has proved to be more 
difficult as explained below. 

In the presence of turbulence, the magnetic and velocity fields can be decomposed into 
a mean and fluctuating parts: B = (B) + b and V = (V) + v, where the (•) stands for an 
average on the realization of the small-scale fields. Substitution of this decomposition into 
Eq. (CQ) and averaging yield the following equation for the mean magnetic field: 

d t (B) + (V) ■ V(B) = (B) ■ V (V) + r/V 2 (B) + V x e , (2) 

where e = (v x b) is the electromotive force. In the framework of mean-field dynamo theory 
[2] e is linear in the mean magnetic field with the following expansion: 

ei = anjiBj) + f3 ijk V 3 (B k ) + . . . . (3) 

where a and (3 are tensors. The symmetric part of a acts as a source term in Eq. (jSJ), creating 
poloidal (resp. toroidal) field from toroidal (resp. poloidal) one. The antisymmetric part of 
the a tensor is interpreted as a mean pumping. The symmetric part of the a effect has been 
thought to be generated by a helical turbulence, which is likely to be induced by Coriolis 
force in stellar convection zones. This type of dynamo is thus classified as aVt if the Q effect 
is stronger than this effect, or a 2 type if this effect dominates over the Q. In convecting 
stars without pronounced differential rotation, the latter is the only possible mechanism as 



there is no effect. The antisymmetric part of (3 is just a turbulent diffusivity which adds 
up to the molecular one whereas the symmetric part may contain additional source terms. 

When the large-scale magnetic field is sufficiently strong, it begins to influence the ve- 
locity field. At some point it will severely quench the generation coefficient (such as the 
a effect), thereby saturating dynamo action. The effect of an uniform magnetic field on 
magnetohydrodynamical (MHD) turbulence was addressed by [6] who showed that the a 
effect is reduced both for weak and strong magnetic field. However, the effect of fluctuating 
magnetic fields on the a effect remains controversial. It has been argued that the dynamo 



re turbulent transport 
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saturates when (B) oc -R~"a/ (v 2 ) due to small-scale magnetic fields. Here, R m is the mag- 
netic Reynolds number. If the small-scale magnetic field grows preferentially to equipartition 
j3], the large-scale magnetic field must be far below equipartition with the previous formula 
holding with n — 1. Alternatively, if the growth of the small-scale field is limited the 
large scale field can almost reach equipartition and the previous formula holds for n = 0. 

The purpose of this Letter is to investigate a novel saturation mechanism which has 
received little attention: the effect of a stable large-scale shear flow. On the one hand, 
strong shear is good for dynamo as it creates magnetic energy via the Q effect. But, a 
strong shear can reduce turbulent transport via shear stabilization [9]. Indeed, the shear 
has been shown to significantly reduce the turbulence intensity and t 
of angular momentum, particle mixing and magnetic diffusion jlO . 
might have a similar effect on the generation of magnetic field, it is crucial to compute 
self-consistently the a effect incorporating the effect of shear (or differential rotation). This 
is especially the case for solar dynamo which is often envisioned to take place at the base of 
the convection zone, where the shear is quite strong, e.g. to compensate the weakness of the 
interface dynamo 13j. Should the a effect be significantly quenched due to shear, it would 
put severe constraints on the magnetic fields that can be created by dynamo action. For 
instance, the dynamo number which characterizes the efficiency of the afl dynamo is given 
by D = aQL 3 /(r] + f3) 2 {14J], where Q is the differential rotation and L is a characteristic 
scale of the system. In this Letter, we show that the intensity of the a effect depends on 
the strength of the differential rotation, and consequently that the dynamo number is not 
simply proportional to the shear intensity. We also provide a consistent theory of momentum 
transport in sheared magnetized plasmas, necessary to understand large-scale shear flows 
(e.g. radial differential rotation). 
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We consider 3D MHD forced turbulence in an incompressible conducting fluid governed 

by: 



ftV + V'VV = - Vp + B VB + i/AV + f , 

V ■ V = , (4) 

together with Eq ([T]). Here, B is the magnetic field given in the unit of Alfven speed, p is the 
total (hydrodynamical + magnetic) pressure, and f is the small-scale forcing. To study the 
effect of shear flows and magnetic fields on small-scale turbulence, we prescribe a large-scale 
flow of the form (V) = — xQy and a constant large-scale magnetic field (B) = B y parallel 
with (V). Then, to solve the equations for the fluctuating velocity field, v = V — (V), 
and magnetic field, b = B — (B), we use the quasi-linear approximation assuming that the 
interaction between fluctuating fields is negligible compared to the interaction between large 



and smal 
[e.g. see 



-scale fields. Note that this is well justified for weak turbulence due to shear flow 



1(3 ]. As a general solution to small-scale equations is not analytically tractable, 
we consider the two limits of weak magnetic field, where the effect of the shear dominates 
that of the magnetic field, and of strong magnetic field, where the effect of the magnetic 
field dominates that of the shear. Specifically, the case 7 = Bok y /Q <C 1 (weak magnetic 
field) and 7 ^> 1 (strong magnetic field) are considered in the following, where k~ l is a 
characteristic scale of turbulent motions. 

Our main interest in this Letter are the Reynolds stress and the electromotive force, which 
determine the growth/decay of the large-scale velocity field and the large-scale magnetic field, 
respectively. First, the Reynolds stress gives a turbulent viscosity (u x u y ) = v T Q, which adds 
up to the molecular viscosity. Second, for an uniform magnetic field, the electromotive force 
([3]) reduces to the term proportional to a only. Furthermore, as the large-scale magnetic 
field is in the y direction, the electromotive force reduces to e« = a, iy B y . The source term 
for the large-scale magnetic field is given by the curl of the electromotive force, for instance, 
the source term for (B z ) is d x e y — d y e x . Assuming that the variation in the x-direction is 
much larger than in y as is the case for the large-scale flow (V)(x) (for instance in the thin 
solar tachocline, the gradient in the radial direction is the dominant one), the main source 
of poloidal magnetic field is d x e y . Consequently, we focus here on the y component of the 
electromotive force 

£ y = (v z b x - v x b z ) , (5) 



and, in the following, we set a = a yy and refers to it as the alpha effect. 

To calculate the correlation functions involved in the transports coefficients [see Eq. (jSJ)], 
we consider an incompressible forcing which is spatially homogeneous and temporally short 
correlated with the correlation time r/. Specifically, in Fourier space, the correlation of the 
forcing is given by: 

(fi^tjf^h)) = T f (2n) 3 5(k 1 + k 2 )5(t 1 -t 2 )x 

<Mk 2 ) , (6) 

where the tilde denotes a Fourier-transform with respect to the spatial variable. As noted 
previously, the a effect can be linked to the helicity of the turbulent flow. Consequently, we 
consider a forcing with both a non-helical part (with energy spectrum E) and a helical part 
(with helicity spectrum H) given by: 



(k) = E{k) (5 lm - + ie lmp k p H(k) . (7) 



It is important to note that in the absence of shear or magnetic field, the kinetic energy and 
helicity of the homogeneous flow v driven by the forcing become: 

,2 



e 



4r f 

h = (v ■ V x v ) - 



/ dk—H{k) . 



(2vr) 2 

Weak magnetic field - To solve the equations for the small-scale fields in the weak magnetic 
field regime, we expand any field ip in powers of 7 = B k y /Q, = ip + 7^1 + • • • an d 
solve order by order. The leading term and the first correction to the turbulent transport 
coefficient are obtained from the first four terms in the preceding expansion. The turbulent 
viscosity is found to be even in Bq whereas the electromotive force is odd: v T = + + . . . 

and a — a\ + 013 + 

To investigate the effect of shear we focus on the limit of strong shear where the shear is 
stronger than the diffusion, characterized by the small parameter £ = uky/Q -C 1. Note that 
this limit is relevant in astrophysical applications such as in the solar tachocline where the 
shearing rate is much larger than the diffusion rate for reasonable values of the parameters. 
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In this limit (£ < 1), the turbulent viscosity can be obtained, after a long algebra, as: 



k 2 k 2 

"F + | T < k > 



Here, k 2 H = k 2 + ^ and T(k) = ^1^1/2^ — arct an (/c^ //;;#). Eq. ([9]) is the kinematic 
result [Bq = 0) and shows that v T is strongly reduced by the shear with scaling Q~ 2 (see 
also [10]). In the 2D case (k z = 0), v T is negative (inverse energy cascade), whereas, in 
3D the second term in Eq. (Q dominates over the first one, making the turbulent viscosity 
positive. It is interesting to note that v T is proportional only to the energy part of the 
forcing, independent of the helical part. The first correction term due to the magnetic field 
to the kinematic result is obtained at second order in 7 as follows: 

T/jB ° 2 f d 3 kk 2 H E(k)T(2/3) ( ^ x 



9(2?r) 3 fi 4 J n \2U 
T(k) 2 



3*3 



Here, T is the Gamma function. The correction (flOl) scales as B 2 tt~ 10 / 3 and can be positive 
or negative depending on the values of the parameters. In 2D (k z = 0), Eq. (1101) is obviously 
positive, slowing down the inverse cascade. In comparison, in 3D, the second term in (flOl) 
dominates the first making v\ negative. In both cases, the correction term (fTUj) is always of 
the opposite sign to the leading order term ([9]). Comparing Eq. fflOl) and Eq. ([9]), we find a 
crossover scale L y at which z/J in Eq. (flQl) becomes larger than : L y Q^BqU^ 1 . 

Similarly, the leading order contribution to the electromotive force is found to be propor- 
tional to B and, consequently, the a effect to leading order is independent of B . It is the 
kinematic a effect which can be computed by ignoring the back-reaction of the magnetic 
field on the flow. In the strong-shear limit 1) this gives: 

'1/3 



ai ~ - whI d3kk2H ^ r( -¥— 



vk l 

Therefore, to leading order, the a effect is reduced by shear proportionally to f2~ 5//3 . This 
result shows, for the first time, that the a effect can significantly be reduced by a strong 
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shear. In contrast to the turbulent viscosity, the a effect is proportional only to the non- 
reflectionally symmetric part of the forcing. This agrees with the expectation that the a 
effect is present only for flow with helicity, which results from the helical forcing with helicity 
spectrum H in our case. 

To next order, we compute the correction to the a effect due to magnetic field. In the 
strong shear limit C 1) we can then obtain the following: 

J * k ** H i k ™h~*&- (12) 

This correction has obviously the opposite sign to the leading order term (ITT1) . Therefore, 
both the magnetic field and the shear quench the a effect. The ratio of these two terms 
is of the same order of magnitude as the ratio of the first correction to the leading order 
term in the turbulent viscosity. In other words, the a effect is suppressed by magnetic fields 
when Bq ~ Q/kRll 3 , where R m = Q/r]k 2 is the magnetic Reynolds number. This is an 
important result showing that a sufficiently strong shear can modify the efficiency of the 
dynamo saturation, leading to a completely new scaling. This a quenching due to flow shear 
has not been investigated in previous works which did not incorporate the effect of shear 
or incorporated it only perturbatively suggesting the critical magnetic held to be either 
independent of R m [8[ or to scale as R^ |7j]. 

Strong magnetic field.- In the case of a strong magnetic field (7 ^> 1), we obtain a WKB 
solution for the small-scale fields and then compute the turbulent transport coefficients with 
the help of Eq. ([6]). First, the turbulent viscosity is given by: 

Jd s kE(k)^x (13) 



T T f 

V - 



(27T)' 



kl k 2 z kyBl 



4B 2 k 2 y k 2 k 2 [B 2 k 2 + v 2 k A Y 

to leading order in 7. Assuming isotropic forcing and performing the integration on the 
angular variables, we find that, for strong magnetic field (B* ^> 1, where Z?* = iBol/vk), the 
turbulent viscosity is reduced by B as B^ 1 due to alfvenization of turbulence 12|. Note 
that in 2D (k z = 0) v T is positive, in agreement with [12 ; j. Furthermore, the second term in 
Eq. (TIBl) is also positive. Therefore, v T in 3D MHD is positive with larger value compared 
to the 2D MHD case. 

Similarly, the electromotive force can be computed to leading order, with the following 
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result for the a effect: 



The integration of Eq. ( TI41) over the angular variables gives: 

arctan(.B* 



"0 = 777 v7 / dk 



(2VT)2 J Q V 2 B 2 



1 - 



Therefore, for strong magnetic field (B* ^> 1) the a effect scales as B~ 2 . For a forcing with 
finite correlation time a different scaling, B~ n with 2 < n < 3, can be expected 0, fig . 

While the formula (jT3j) shows that in the strong magnetic regime the main quenching of 
the a effect is due to the magnetic field, shear flow is found to have an interesting effect on 
a at second order in 7 _1 as follows: 

a 2 = ^l— [ dk ^§ (-^] x (16) 
2(2n) 2 J v 2 Bt \vk 2 ) 1 ; 



2 arctan^*) 



For strong magnetic field (5* 3> 1) we see in ([TBI) that the first correction is of the opposite 
sign to the leading order and that it scales as Q 2 B~ 4 . Comparing (TIB"]) with (|T5|) we can easily 
see that the correction (IT51) is always negligible compared to the leading order term ( fl6l) as 
B k/Q ^> 1 in the strong magnetic field regime. Consequently, the a effect is suppressed 
when £>* 3> 1, i.e. for £> 3> Vt/kR m . 

Discussion - We show an a-quenching due to shear (a oc for the first time. Depending 
on the ratio of the large-scale magnetic field to the shear, we found that the a effect is 
quenched by large-scale magnetic field with different scalings with R m . Specifically, the 
critical magnetic field strength B c , above which the a effect is suppressed by magnetic field, 
is B c ~ Q/kRlr( 3 for 7 = B^k/Vt ^> 1 while B c ~ Q/kR m for B$k <C Q, where k~ x is the 
typical small scale of turbulence. For instance, in the solar tachocline (fi ~ 3 x 10 _6 s _1 , 
v ~ 10 _2 m 2 -s _1 and B ~ 6m-s _1 ), the strong magnetic field regime (B Q k Q) is valid on 
scales less than 2 x 10 6 m, which is approximately the size of the tachocline. Therefore, not 
only magnetic field but also shear can dramatically quench the a effect hindering large-scale 
dynamos. 

The turbulent viscosity is reduced by strong magnetic field (with a scaling B~ 2 ) while 
in the opposite limit of weak magnetic field is quenched by shear (with a scaling Q~ 2 ). In 
both cases, the turbulent viscosity is positive with a larger value in 3D than in 2D. In the 
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weak magnetic field case, magnetic field becomes important on scale L y < Q^BqU^ 1 . In 
case of the Sun, L y < 5.4 x 10 15 m, suggesting that turbulent viscosity is severely quenched 
by magnetic field. 

These results have crucial implications for dynamics and angular momentum transport 
in astrophysical or laboratory plasmas which are often envisioned to be efficient. In partic- 
ular, quenching by shear and/or magnetic field should be incorporated when assessing the 
efficiency of dynamo, e.g. the dynamo number D. For instance, our result a oc f2~ 5 / 3 makes 
D oc f2~ 2 / 3 , decreasing with strong shear rather than increasing (D oc Q) as previously 
thought. Furthermore, the nonlinear dependence of turbulent viscosity on Bq can offer an 
interesting mechanism for the time variability such as torsional oscillation in the Sun 16]. 
It will be interesting to extend our theory to incorporate the effects of rotation which will 
consistently give rise to a effect and non-diffusive momentum transport, (A effect), due to 
shear-induced anisotropy [l7]. How the A effect, a effect, turbulent viscosity and particle 
transport are affected by rotation, magnetic field and shear would be of great interest with 
important implications. The controversial issue on the quenching of j3 effect in 3D MHD 
should also be investigated by incorporating the gradient of large-scale magnetic field. These 
issues will be addressed in future publications. 
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